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TEST PARTICLE PROPAGATION IN MAGNETOSTATIC TURBULENCE 


n, THE LOCAL APPROXIMATION METHOD 


ABSTRACT 

A new approximation method for statistical mechanics is presented and 
applied to a class of problems which contains the test particle propagation 
problem of paper I. All of the available basic equations used in statistical 
mechanics are cast in the form of a single equation which is integrodifferential 
in time and which is then used as the starting point for the construction of the 
"local approximation method". Simplification of the integrodifferential equation 
is achieved through approximation to the Laplace transform of its kernel. 

The approximation is valid near the origin in the Laplace space and is based 
on the assumption of small Laplace variable. No other small parameter is 
necessary for the construction of this approximation method. The n'th level of 
approximation is constructed formally and the first five levels of approximation 
are calculated explicitly. It is shown that each level of approximation is governed 
by an inhomogeneous partial differential equation in time with time independent 
operator coefficients. The order In time of these partial differential equations 
is found to increase as n does. At n = 0 the "most local" first order partial 
differential equation which governs the Markovian limit is regained. As n 
reaches infinity the original integrodifferential equation which is non-local in 
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time is recovered* The gap between the Markovian limit and the non-local 
integro differential equation is found spanned, with increasing n, by a system 
of increasingly non-local theories. The local approximation method is applied 
to two examples. The first of these is a simple exactly soluble case which is 
obtained through the choice of an appropriate C-number kernel in the integro- 
differential equation. The exact solution and all of its local approximations are 
compared for this case. The second example is the quasi-linear diabatic 
equation of paper I. which governs the probability distribution function for test 
particles in a simple model of magnetostatic turbulence. In this case the first 
(n = 1) local approximation is evaluated numerically and presented for the same 
initial conditions which were used in paper I. For both examples, with the ex- 
ception of the Markovian limit, the local approximations generally produce very 
accurate approximations to the exact solutions whenever they evolve slowly in 
time, or they give time averages of the exact solutions whenever they evolve 
rapidly. In particular, the failure demonstrated in paper I. of the Fokker- 
Planck {velocity space diffusion) equation which governs the test particle prob- 
ability distribution function in the Markovian limit is removed by the "almost 
Markovian" first local approximation. 
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TEST PARTICLE PROPAGATION IN MAGNETOSTATIC TURBULENCE 
II. THE LOCAL APPROXIMATION METHOD 

I. INTRODUCTION 

In the first paper in this series^ we have seen that the familiar Markovian 
characterization of the propagation of test particles in a simple model of mag- 
netostatic turbulence fails to accurately represent the quasi-linear approxima- 
tion to the probability distribution function. In our effort to replace the adiabatic 
approximation, which leads to the Markovian picture, we have discovered a new 
’local approximation method”. In this paper we develop the local approximation 
method as it applies to a class of integrodifferentiai equations which contains 
the quasi-linear diabatic equation that we constructed in paper I. as well as all 
of the fundamental equations which govern other particle and fluid systems at 
the kinetic and transport levels. 

We outline the local approximation method in general terms, construct the 
n’th level of approximation formally, and calculate the first five levels of ap- 
proximation explicitly. We show that each level of approximation is governed 
by an inhomogeneous partial differential equation in time with time independent 
operator coefficients. At n = 0 the adiabatic approximation is recovered. 
Because the adiabatic approximation is governed by a first order partial differ- 
ential equation in time, its evolution at any instant depends only on its state at 
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that same instant. This first order equation gives the differential analogue of 
the Markov chain in which the state of a physical system is imagined to proceed 
in distinct steps according to a transition probability which is independent of 
time. In the Markov chain the state of the system after a step forward in time 
can be given entirely in terms of the state of the system at the preceding step. 

In the Markovian picture, the adiabatic approximation, the solution of the 
system under consideration is governed by 'local" equations; the evolutioa of 
the system at any instant does not depend on the state of the system at previous 
times. 

As n is increased the order in time of the partial differential equation which 
generates the n'th local approximation also increases. At n = co, the local ap- 
proximation is governed by the infinite order partial differential equation which 
is equivalent to the original integrodifferential equation; i.e., the exact solution 
is regained as n goes to infinity. This integrodifferential equation is non-local 
in time; the evolution of its solution at any instant depends on the state of its 
solution for all previous times. 

Within the local approximation method, the zero'th local (or adiabatic) ap- 
proximation produces the "most local" Markovian description of the solution 
under consideration. The other extreme (n = <^) gives the exact solution which 
is governed by the integrodifferential equation which is non-local in time. 
Between these two extremes we have an infinite set of partial differential equa- 
tions in time with members whose order increases with increasing n, and 
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therefore, whose degree of "localization" decreases with increasing n. The 
local approximation method provides a mechanism for bridging the large gap 
between the completely Markovian theory and the completely non-Markovian 
theory with a system of increasingly non-local theories. 

The construction of the local approximation method requires no small 
parameter. Consequently, we feel that this method will provide many important 
generalizations of presently available kinetic theories into parameter regimes 
which have not been theoretically accessible previously. Simplifications of the 
original intergro differential equation are produced by constructing its Laplace 
transform and then approximating the Laplace transform of its kernel near 
the origin; i.e., it is the small Laplace variable which is used to produce the 
approximation. Fade approximants^ to the Laplace transform of the kernel 'are 
used to insure that the approximation which is constructed for small Laplace 
variable does not fail catastrophically for large values of the Laplace variable. 
Our application of this method is limited in this paper to those approximants 
which lie near the diagonal of the Fade table and are equivalent to continued 
fraction convergents.^ The order of the convergent for the kernel transform 
determines the level of the local approximation to the solution. We show that 
within the class of integrodifferential equations which we consider, there is an 
identifiable sub-class which contains members equivalent to finite order partial 
differential equations in time with time independent operator coefficients. This 
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sub-class is identified by those kernel transforms which can be expressed as 
the operator ratio of two finite order polynomials in the Laplace variable.^ 

In this sub“class the local approximation method converges to the exact solution 
for finite n. All higher order local approximations are equivalent, and produce 
the exact solution. 

The source terms in the inhomogeneous equations which generate the local 
approximations are zero for all time except at the origin in time where they 
are singular. Thus, the solutions of these equations are double valued at the 
origin in time; they jump from the true initial conditions to "reinitialized" 
conditions. Each local approximation has its own unique reinitialization. At 
least in the examples which we will present in this paper, we often find rapid 
transients in the exact solutions which are short lived and which lead to slow 
evolution for larger times. The low order local approximations which we 
investigate explicitly in this paper are not capable of supporting these rapid 
transients. Nevertheless, reinitialization seems to choose the appropriate 
particular solution so that once the transients are over the local approximations 
and the exact solution agree very well. During the transients the local approxi- 
mations give an indication of the time average behavior of the exact solution. 
Thus, with the local approximation method we are able to study the initial value 
problem without the necessity of computing the details of the often intricate 
initial transients. 
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We apply the local approximation method to specific problems and present 
the results. In section IH, we consider an exactly soluble example of the sub- 
class of integrodifferential eq.uations discussed above. A simple kernel function 
is chosen which leads to the reduction of the original integrodifferential equation 
to an ordinary third order differential equation with constant C -number coeffi- 
cients. Thus, the exact solution IS readily available. The system of local 
approximations truncates at n = 4 in this ease, and produces a series of first, 
second and one third order (the exact) ordinary differential equations with 
constant C-number coefficients. Thus, each of the local approximations is 
readily available as well. The various local approximations are presented and 
compared to the exact solution. The details of this comparison are presented in 
section IE. Generally spealting, we find that when the local approximation 
method produces meaningful results they either take the form of very accurate 
approximations to the exact solution whenever it evolves slowly in time, or they 
give a time-averaged representation of the exact solution if it evolves rapidly. 

In section IV, we apply the local approximation method to the test particle 
propagation problem. With our choice of exponential correlation function given 
in paper I., the system of local approximations to the quasi-linear diabatic 
probability distribution function truncates at n = 4. In paper I. we have presented 
several numerical evaluations of the exact quasi-linear diabatic solutions as 
well as the adiabatic approximation to one of these solutions for an important 
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initial condition. In the terminology of this paper we have presented several 
numerical evaluations of the fourth (exact) and the zero'th (adiabatic) local 
approximations to the quasi-linear diabatic probability distribution function in 
paper L. In this paper we present numerical evaluations of the first local ap- 
proximation to the quasi-linear diabatic solutions for the same initial conditions 
which were considered in paper I. 

In terms of numerical integration, we find that the local approximation 
represents a considerable simplification. Various numerical instabilities which 
appear in the quasi-linear diabatic numerical evaluations never appear in the 

2 3 

approximate numerical evaluations, and a factor of 10 to 10 less computer 
time is found necessary to achieve the same point in time in the approximate 
solutions compared to the exact ones. 

We show by direct comparison that following the initial transients in the 
exact solutions, which last somewhat less than one Larmor period, the first 
local approximation becomes very accurate for each of the initial conditions 
considered. This accuracy persists through the ends of each of the numerical 
evaluations. Some indication of the accuracy of the first local approximation 
for time beyond the end of our numerical evaluations is obtained by a study of 
the Laplace> model expansion of the probability distribution function which was 
introduced in paper I. We show that the first local approximation to some of 
the low order Laplace modes is very accurate for small values of the Laplace 
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variable. We also show that the eigenvalues are accurately given by the first 
local approximation. These results indicate that the accuracy of the first local 
approximation persists. In the next paper in this series we develop this reason- 
ing much further in order to obtain, strong evidence that the first local approxi- 
mation remains accurate for all post transient times. 


II. CONSTRUCTION OF THE LOCAL APPROXIMATION METHOD 

We present the local approjdmation method through its application to the 
equation, 


- e 
'dr 


f d^K(A)f(r-M 


(E.1) 


*'0 

in which the kernel, K(\), is a linear operator which we will assume has a 
Laplace transform with a Maclaurin series expansion in the Laplace variable. 
Other variables (e.g., phase space variables) will be assumed in both K(^) and 
f(r) but will not be exhibited explicitly here. The quasi-linear diabatic equation 
for test particle propagation in magnetostatic turbulence which was constructed 
in paper I. is of this form, and in addition, all of the basic equations used in the 
field of statistical mechanics and some important macro-turbulence theories 


fall in this classification as well (see Appendix A). 

In equation H.l, e is a parameter which is ordinarily assumed small in 
order to produce kinetic approximations for the long time, near equilibrium 
evolution of f(r). ’ The necessity of this assumption normally limits the avail- 
able kinetic theories to low density and/or weak interaction systems for which 
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the time scale which characterizes the duration of a particle collision is much 
smaller than the time scale over which f(r) evolves. In many situations, even 
though these time scales may be well separated for typical particle parameters 
they may become mixed for some more or less important regions in the available 
particle phase space. In paper I. we have seen an example of mixing of the 
interaction and evolution time scales in an important region of phase space 
(near^t = 0) and the consequent failure of the standard adiabatic approximation. 

In this paper we will see that the local approximation method does not require 
the clear cut separation of time scales and is thus able to remove the non- 
uniformity found in paper I. Of further significance, we will find that the con- 
struction of the local approximation method does not require e small. Con- 
sequently, we feel that the local approximation method will provide many im- 
portant generalizations of presently available kinetic theories into parameter 
regimes which have not been theoretically accessible previously. 

Equation II. 1 can be put in the renewal form. 


r-r 


f(r)= f(0) + e 


dX Jy(X) f (r - \) 


^0 


II.2 


in which. 


D(\) = 


ds K(s) 


*^0 


II.3 


8 



E K(\) is a C-number kernel, then equation II. 2 is a Volterra integral equation 
on the semi -infinite domain with a specialized form for the kernel (convolution 
integral under Laplace transformation). Of course, if we are faced with the 
solution of an equation of this form, we can obtain equation II. 1 through differ- 
entiation with respect to time whether or not the kernel is a C-number kernel. 

Our local approximation method, as we have developed it so far, consists 
of the following steps; First we put the equation, or system of equations, under 
consideration into the form of equation n.l. Second, we Laplace transform that 
equation. This Laplace transformation is possible because the integral term 
in equation n.l is a convolution integral under the Laplace transformation 
(however, see Appendix A where the convolution form is dropped). Third, we 
construct a Maclaurin series expansion of the Laplace transform of the kernel 
in powers of the Laplace variable. A basic requirement which we have had to 
assume so far in our construction of this method is that the Laplace transform 
ot the kernel has a Maclaurin series expansion. Fourth, from the Maclaurin 
series expansion we construct an operator continued fraction expansion of the 
kernel. Fifth, we truncate this continued fraction expansion at the n'th order 
in order to produce the n'th "convergents and then we substitute the n'th 
convergent back into the Laplace transform of the original equation in place of 
the exact Laplace transform of the kernel. At this point we have constructed 
the equation in Laplace space which generates the n'th local approximation 
to the Laplace transform of the exact solution. Finally, we inverse Laplace 
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transform this equation in order to obtain the equation which generates the n'th 
local approximation to the exact solution of equation II. 1 in r . 

At the end of the procedure outlined above we have constructed a finite 
order, inhomogeneous, partial differential equation in time with operator co- 
efficients which are constant in time. Except in some special situations which 
will be discussed below, the order of this partial differential equation can range 
from first, when n = 0, to infinite (the integrodifferentiai equation vdiich generates 
the exact solution) when n reaches infinity. The zero'th local approximation, 
which is generated by a first order partial differential equation in time, is 
identical to the adiabatic approximation which was discussed in paper I. Thus, 
the local approximation method produces a series of equations which generate, 
in the lowest order, the most local, Markovian, approximation to the exact 
solution, and in the highest order, the exact solution which is generated by the 
non-local integrodifferentiai equation; the finite orders in between bridge the 
gap between these two extremes. By construction, all of these equations 
generate solutions vdiich have Laplace transforms that accurately approximate 
the Laplace transform of the exact solution when the Laplace variable is small. 
We will see in this paper through the various examples which will be presented, 
that this feature of the local approximations usually results in accurate approx- 
imations to the exact solution whenever it evolves slowly, or, whenever it evolves 
rapidly, to time averages of the exact solution. 
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In the remainder of this section we present a detailed construction of the 
local approximation method as outlined above. 


The Laplace transform of equation 11,1 is, 


pf(p) - f(0) = eK(p) f(p) 


and, in the Laplace space, its exact solution can be written formally as, 


f (p) = [p - eK(p)] ^ f (0) 


The kernel is assumed to have a Maclaurin series expansion in the Laplace 
variable; 




in which the M are the n'th moments of the kernel, defined through. 


= dX. X" K(X) 


For the remainder of this paper we adopt the operator ratio convention, 

Y 

f\j 

With this convention a continued fraction expansion of K(p) can be constructed 
from the Maclaurin series expansion. 


K(p) = + ■ 


1 + ••• 
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in which some of the low order coefficients are. 


ao = Mo 

a2=lM“iMo ^3 -1 m 7^M2 n.lO 

^ - 2 ' [i^ “^“4 - ^ (M?l*l3)=] 

We obtain the n'th convergent, K„(p), to K (p) when we truncate the continued 
fraction expansion in equation II.9 beyond a^. We obtain the n'th local approxi- 
mation to f(p), which we denote, f^ (p), by substituting (p) in equation n.5 for 
K(p). Thus, the equation which generates the n'th local approximation to the 
Laplace transform of the probability distribution function is, 


pT (p)- f(O).. eK>) f^CP) 

Any kernel transform which, can be expressed as a ratio of an I'th order 
to an m'th order polynomial in p can then be expanded as in equation II.9 with 
all a^ = 0 for n > 1 + m + 1. All convergents of higher order than 1 -f- m + 1 
are identical, and are exactly K{p). Thus, if the Laplace transform of the kernel 
can be expressed as the operator ratio of two finite order polynomials in p, 
then the local approximation method must converge on the exact solution for 
finite n. 
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The (L/M) Pade^ approximant^ to K(p) is given by, 

(L/M')=Jl 11.12 

Qm(p) 

m which (p) is a polynomial in p of degree at most L, and (p) is a polynomial 
in p of degree at most M, and in which the coefficients in and Q have been 
chosen so that a power series expansion in p of equation H.12 leads to the first 
L + M + 1 terms in the Maclaurin series expansion of K(p). A Fade approximant 
exists for each combination of L and M. Thus, an array can be constructed, 
called the Fade Table, in which the element of the array in the L’th row and 
M'th column is given by (L/M). A portion of this array is given in figure 1 for 
the special example in which K(p) is the ratio of a first to a second order poly- 
nomial in p. In this case, all elements of the Fade" Table with L > 1 and M ^ 2 
are identical and equal to K{p). The M - 0 column contains the L*th order 
truncations of the Maclaurin series expansion of K(p); some other elements of 
the table are exhibited. The boxes with dark borders contain the n'th convergent 
to K(p) for n = L + M . For the special example being considered here, all 
convergents with n > 4 lie in the infinite block of identical approximants given 
by K(p). The n'th local approximation to f(p) is obtained by substituting the 
(L/M) approximant which has L + M = n, and which lies in the "staircase" of 
dark boxes in the Fade' Table, for K(p) in equation II.5. 

We have not investigated the consequences of choosing approximants to K(p) 
other than those which correspond to the convergents. However, we point out 
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Q-NUMBER PADE TABLE 



I + iM2"'M3 p 

Figure 1. The Fade table for a Q-number kernel whose Laplace transform can be expressed 
as the ratio of a first to a second order polynomial in the Laplace variable. The operator 
ratios shown here are to be interpreted according to equation II. 8 in the text. 







the connection between the convergents and the Fade approximants to indicate 
a possible generalization of the particular construction which will be presented 
here. Fade approximants have been used to accomplish approximate Laplace 
transform inversion, and they have also played an important role as an approxi- 

4 

matron technique in Hilbert space for various quantum-mechanical models. 

The equation which generates the n'th local approximation to f(r) can -be 
constructed by inverse Laplace transforming equation 11,11. By inverting this 
equation directly, we find, 

r 

\ (^) f„(T - >.) IL13 

in which the (X) are the inverse Laplace transforms of the convergents for 
the kernel transform. In Appendix C some of the are computed explicitly. 

Here we will show that the fn("^) also obey equivalent inhomogeneous partial 
differential equations in time with time independent operator coefficients. These 
equations are obtained by first manipulating equation n. 11 into a convenient form 
and then performing the inverse Laplace transform. This procedure is outlined 
in the following for n = 0 through 4. 
a), n = 0 

For n = 0, we have, 



Pfo(p) - f (0) = 


11.14 
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and therefore. 


^^0 
3r ■ 


11.15 


with the proviso that fo(0) — f(0)> The formal solution of equation 11,15 is. 


= f(0) 


n.i6 


Equation 11.15 can be compared with equations 1.26 and 1.27 or equation 1.36 of 
paper I. in order to see that the zero'th local approximation (n = 0) is just the 
adiabatic approximation. In this approximation we replace the original integro- 
differential equation (in time) with a first order differential equation (in time). 

We have seen in paper I in our study of test particle propagation in magnetostatic 
turbulence that this level of approximation produces the Markovian description 
of f (t) which is also available from a Fokker-Planck equation in velocity space 
for f(r). In analogy with the Markov chain, we feel that this description of f ("^i) 
is as local a description of its evolution as can be constructed. We use ”local” 

In contrast to equation II. 1 which contains an obvious memory of the past and is 
clearly non-local in time, 
b), n = 1. 

The first local approximation is generated by, 

pf,(p) - f (0) = 6K,(p) f/p) = e(Mo - pM,) f, (p) IL17 


16 



The inverse Laplace transform of this equation is. 


3f 


[1 + eU^] = - eS(r) M^f(O) 


11.18 


which has a formal solution given by, 


flC-r) = e ^ [1 + €U^\ i 


( 0 ) 


11.19 


Notice that, 

f (r) ►[! + eM ] ^ f (0) H.20 

Because of the presence of the Dirac delta function on the right hand side of 
equation 11.18, f^lr) is discontinuous at r = o, jumping from the initial condition, 
f(0), to the "reinitialized" condition given by equation 11.20. Equation 11.18 is 
also a first order differential equation in time, and in this respect, is similar 
to equation 11.15. However, the presence of the Dirac delta function enforces a 
memory of the origin in time in equation 11.18 thereby leading to a non-Mar kovian 
kinetic theory for f^(r). In a sense, ‘this level of approximation is local com- 
pared to the integrodifferential equation which generates the exact solution, 
but not as local as the Markovian description, 
c). n = 2, 

The second local approximation is generated by, 

P^2 (P) " f (0) = e KjCp) ?2(p) 11.21 
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In order to construct the equation which generates 


K 2 (p) = ap + 3j^p(l + ajP)"^ IL22 

must be first rewritten in an appropriate way. 

KjCp) = Bq + [(1 + 32P) P 

= ag -s- 4 p)] ^ p jj 23 

= ao + (p + a^a"^p 


By substituting this expression into equation II.21 and operating on the equation 
(P 3^1 3-2^ 3-1^ )) obtain, 


(p + a^a'ia^i) [p'f2(p) " f(0)J = ->■ ^2(P) 

This form can easily be inverse Laplace transformed to 

^+^2 ^ + A2f2 = eS(r)D2f(0) 

in which, 

^2 = - 

and 

Dj = 3(5 + a^aj^ 


11.24 

( 


11.25 


11.26 

11.27 


IL28 


18 


&BlPEOpUCIBILirY OF THE 
©R5MNAL PAGE IS POOP 



In this case reinitialization leads to 


f2(0-^)= f(0) 


"Sr 




11.29 


Higher order local approximations can be obtained in a manner similar to 
that used above for the second local approximation. As higher order convergents 
are considered the operation given by equation 11.23 must be repeated an increas- 
ing number of times in order, in each case, to reach an equation for (p) which 
is similar to equation 11.24. In this equation, only powers of p multiply either 

r\j 

f^ (p) or f(0). Once this form is reached the inverse Laplace transformation can 
be accomplished easily, 
d). n = 3. 

We find that the third local approximation is generated by, 

C3 ^ + B3^ + A3f3 = eS(r)D3f(0) + eS(T) E3f(0) n .30 

with, 


^3 = - 

ea„a,(l -i- 

IL31 


ai(l + 3332^)"^ 32 ^ 3 ^ - - ^^0 

11.32 


1 _ ea^(l + 3283 ^ I"’- 

11.33 

Da = 

ao + ai(l 4 - 333;!)"^ a’^ 

11.34 
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and, 


-iN-l 


E 3 - a^(l + 83 83 ) 


and with reinitialization giv_en_by., 


f3(0") =c-if(0) 


313 ( 0 -^) 

3t 


= €3^(03-83)8-1 -B3C3I] f(0) 


e). n = 4. 


The fourth local approximation is generated by, 


33 f 32 f. 3 r 

‘i 4 n 4 


+ c, 1 + B, f A, = e S (r) DJ (0) + £ S (r) E, f (0) 


3 t 3 Br2 ^ 


with, 


A4 = £8^8384 33^32 


B4 = a Ja3a;ia“ia;ia7 - £ (a3a4ia'ia;i + a3a;ia3ia-i + a^^a^^) a„ - £a3a;ia;ia 

I 

C4 = 3^338418"! f 9394133^ + a"i] a"i - e [8 q + a^Cags"! + 1) a~i] 

D4 = - e [1 - 83(838-13-1 + 838-18-1 + a^l) a-la^ - a^Sga^la-la'l] 

and 

E 4 =a^ +a^(l + a 3 a;i)a;i 


IL35 

11.36 


11.37 

IL38 
2^3 11.39 
11,40 
n.41 

n.42 
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and with reinitialization given by. 


f^(O^) = f(0) 

— = eD f(0) 11.43 

dr ^ 

B2f (0+) 

= e [K(0) - C^D, - E4] f (0) 

The procedure outlined above could be carried out further in order to 
obtain even higher order local approximations. In the remainder of this paper, 
however, we will apply this method to two examples in which the fourth and 
all higher local approximations generate the exact solutions. For our purposes 
here, it is not necessary to compute the higher order approximations. 

As we proceed down the staircase in the Fade" table with increasing n, or 
equivalently, as we adopt higher order convergents for the kernel transform, 
we produce higher order local approximations for f(T). These local approxi- 
mations are generated by partial differential equations of increasing order in 
time (the order is given by M + 1) with operator coefficients which are constant 
in time, and with reinitializations appropriate for each approximation. If the 
kernel transform can be expressed as the ratio of two finite order polynomials 
in p, then this procedure truncates at a finite order, and in that order it produces 
the partial differential equation which is equivalent to the original integrodif- 
ferential equation. Otherwise, this procedure converges on the exact solution, 
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which is generated by a true integrodifferential equation, only as n goes to 
infinity. 

All of the partial differential equations which generate the various local 
approximations have one other feature in common; they all generate solutions 

I 

which have Laplace transforms that accurately approximate the Laplace trans- 
form of the exact solution when the Laplace variable is small. We view this 
feature as the guiding principle in the construction of the local approximation 
method, and in this sense, consider the particular construction which we have 
presented here as only one of many possible such constructions with the same 
guiding principle. 

We have seen that the zero'th local approximation is identical to the adiabatic 
approximation. At this level of approximation we reproduce the Markovian 
description of the probability distribution function. In general, at this level the 
local approximation is governed by a first order differential equation in time 
which has no reinitialization. This equation is time-translation invariant {it 
does not contain a memory of the origin in time), and it produces changes of 
f^j(r) at any instant in time which depend only on fo(r) at that instant, very much 
in the spirit of the Markov chain. We call this equation a local equation in 
contrast to the non-local integrodifferential equation. The local approximation 
method, therefore, provides a bridge from the most local Markovian description 
to the completely non-local integrodifferential description. In paper I. of this 
series we have seen an example in which the Markovian description fails to 
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produce an accurate approximation to the exact solution. In section IV. of this 
paper we will reexamine that case and show that the "almost Markovian" fj(r) 
provides a very good approximation to the exact solution. 

III. APPLICATION TO AN EXACTLY SOLUBLE PROBLEM WITH A 
C -NUMBER KERNEL 

In the event that the kernel under consideration in equation H.l is a C- 
numher function, the equations which generate the local approximations to f(v) 
reduce to ordinary differential equations with constant coefficients. Thus, their 
solutions are readily available. In this section we present some examples of 
these solutions for the C-number kernel, 

K(A) = cos ^ 11.44 

^ in 'Which is a parameter which lies in the domain, o £ ^ 1. The Laplace 

transform of this kernel is , 


K(P) = —LLf: n.45 

1 4 (P -f 1 ^? 

Since this kernel transform is the ratio of a first to a second order polynomial 
in p, we expect the fourth, and all higher, local approximations to produce the 
exact solution. Indeed, with this kernel, and with the use of repeated time 
differentiation, we can show that equation II. 1 is equivalent to, 

+ 2/x + (1 + _ e) — ~ e^f = 0 jj 46 
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with the additional initial conditions, 


^f(O) - 




11.47 


The kernel transform given in equation n.44 does have a Maclaurin. series 
expansion with some of the lower order moments given by, 


M 


0 


1 



1 

(1 


^ _ 2^(3 - ^ _ 6(1 - 6/j? ■}■ 


11.48 


M - - 10/j? + 5) 

' ' (1 ^ 

Through the substitution of these moments into equation n.lO, we are able to 
show that equations n.37 through n.43, which generate the fourth local approxi- 
mation, reduce exactly to equation IT. 46, which generates the exact solution. In 
this case the entire system of local approximations is comprised of foCT") through 
5 3-^id the Fade table given in figure 1 is directly applicable with the under- 
standing that the moments which appear there are C-numbers. We have chosen 
to study this particular C-number kernel since it has the same functional de- 
pendence on fjL and X as the Q-number kernel which appears in the quasi-linear 
diabatic equation which was constructed in paper I., and whose first local 
approximation will be studied in detail in the next section of this paper. 
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The local approximations. to f(r), includii^ the fourth local approximation, 
(the exact solution), are plotted in figures 2 through 5 for e =1.0 and several 
values of ix. Notice that in the vicimty oi p. = 1.0 the exact solution contains 
a short-lived weak transient which leads quickly into simple exponential decay. 
For these values of p the various local approximations start with individual 
values and slopes subject to their respective reinitializations, but they all 
quickly become essentially equivalent and then accurately represent the exact 
solution once it reaches the post-transient slow exponential decay phase. In 
the vicinity of p = 0 the transients in the exact solution are no lor^er weak nor 
short lived. Instead, the exact solution oscillates about a mean value which 
decays slowly with time, and with an amplitude which also decays slowly with 
time; both decays are exponential but slow. At p = 0, the exact solution oscil- 
lates with constant amplitude about a mean value which does not decay with 
time. In the region of small p -values (t) takes on what appears visually to 

r 

be a time average of the exact solution and essentially duplicates fj(r) 
except for a very short transient. The adiabatic approximation, f^ (r), provides 
a reasonable picture of the average behavior of the exact solution for large 
values of p , but, for small p both f^ (r) and f^ (r) are definitely preferable 
representations of this average behavior. Only f 3 (r) fails as an approximation 
to the exact solution for some values of p and large times. When p < 0.6, 
f ^(r) contains a positive exponential mode which causes this solution to grow 
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Figure 2. The zeroUh (adiabatic) through fourth (exact) local approximations for the exactly 
soluble problem with a C -number kernel. This is a strong coupling ( e = 1) example with 
jj. = 1.0. The zero'th, first and second local approximations are identical. Notice the 
weak, short-lived transients in the exact solution, following which all of the local approxi- 
mations become essentially identical. 




Figure 3. The zero'th (adiabatic) through fourth (exact) local approximations for the exactly 
soluble problem with a C-number kernel. This is a strong coupling (e = 1) example with 
/I = 0.7. The individual reinitializations are readily apparent in this figure. Notice the 
transients in the exact solution are somewhat stronger than in the previous example. 
During the transients the local approximations give a reasonable picture of the average 
behavior of the exact solution. 




Figure 4. The zero'th (adiabatic) through fourth (exact) local approximations for the exactly 
soluble problem with a C-number kernel. This is a strong coupling (e = 1) example with 
IJ. = 0.1. The exact solution is dominated by slowly decaying oscillations. Both the first 
and second local approximations give the average evolution of the exact solution; the 
zero’th local approximation does not appear as useful. The third local approximation 
contains a growing exponential mode which leads to its divergence with increasing time. 




TIME 


Pigure 5, The zero'th (adiabatic) through fourth (exact) local approximations for the exactly 
soluble problem with a C-number kernel. This is a strong coupling (e = 1) example with 
jj. - 0.0. The exact solution oscillates trigonometrically for all time. The first, second, 
and third local approximations are identical and are exactly the time average of the exact 
solution. The zero’th local approximation gives an indication of the average behavior of 
the exact solution, however, its amplitude is shifted upward from the correct value. 



unbounded with increasing time. This exponential growth rate decreases with 
decreasing /i, however, so that as /i decreases i^{r) remains a viable approxi- 
mation for longer times. At /x = 0 this growth rate goes to zero and f3(r) 
becomes identical to (r ) and f2(r). All three of these approximations then 
become equal to the constant time average of the oscillating exact .solution. 

We feel that these results, and those to follow, do indicate the existence of 
a method for generating kinetic approximations which give excellent fits to the 
exact solution whenever it evolves slowly, or else, give an indication of the time 
average behavior of the exact solution it involves rapidly. This method depends 
on the existence of a formal solution in Laplace space which can be simplified 
considerably through a Maclaurin series expansion. This Maclaurin series 
expansion, with its typically finite radius of convergence, is replaced by Pade' 
approximants in order to continue the approximation to the exact solution beyond 
the radius of convergence to p = “ . The accuracy of the approximation for large 
p does not seem to be too signficant since this domain in p seems to affect only 
the initial short-lived transients in the local approximations . Various Fade" 
approximants give various representations of the solution for short times, but 
they are all equivalent for large times when they produce meaningful results at 
all. In applications of the Fade’ approximants in the literature, where the exact 
solution is not available, the validity of the results is typically judged by com- 
puting the predictions of several approximants and comparing them. On this 
basis, in our application here, we would certainly reject the third local 
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approximation since it diverges radically from the zero'th, first, and second 
approximations which all agree with each other reasonably well. 

There is one further aspect of this approximation technique which is partic- 
ularly interesting. In our construction of the various local approximations we 
have relied solely on the smallness of the Laplace variable; the parameter, ^ , 
was not assumed small. To emphasize this point we have plotted the various 
local approximations in figures 2 through 5 for e = 1.0. In paper I. we have 
presented two examples of the more traditional reasonir^ which leads to the 
adiabatic approximation as a leading approximation in an asymptotic expansion 
in the small parameter, e . Here we see that the adiabatic approximation, and 
particularly the higher order local approximations, still jneld information about 
the average behavior of the exact solution even when e is not small. Equation 
II. 1 can be considered a model of the kinetie theory problem for a variety of 

7 8 

physical systems. ’ When e is small we are in the "weak coupling" regime 
for all of these kinetic theories. It appears that with the use of the local approxi- 
mation method we should be able to learn something about the average behavior 
of these other physical systems in the strong coupling r^ime in which the 
particle-particle or particle-wave interaction strength and/or the particle 
density is not assumed as 5 miptotically small. 
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IV. APPLICATION TO THE QUASI-LINEAR DIABATIC EQUATION FOR 
TEST PARTICLE PROPAGATION IN MAGNETOSTATIC TURBULENCE 
In this section we compute the first local approximation to the probability 
distribution function which is the solution of the quasi-linear diabatic equation 
constructed in paper I. for test particle propagation in magnetostatic turbulence. 
As a matter of convenience, we treat equation 1.24 from paper I. which we 
reproduce here: 


= V(i 


f dX.J{0x, X) 


•'o 


r - k) 
'd/j? 


11.49 


in which, 


CM 


N(ja, t) = 


dfM‘ fOLi'. T) 


11.50 


Thus, in this section we obtain numerical solutions of equation 11,18 with e 
replaced by ■ 772 , and with ) given by, 


K(\) = (1 k) 


_!L 


n.5i 


where. 


k) = R(jj.k) cos k 


11.52 
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and where R(^) is the two point correlation function associated with the as- 
sumed homogeneous magnetostatic turbulence. As in paper I. we adopt, 


E(0 = 


11.53 


Then, 


K(p) = (1 - p) 

'htj? 


n.54 


in which, 

P) = 

1 + (P + Ip-I f 

The probability distribution function is obtained by numerically differentiating 
N(ti, r) with respect to jx. 

We have chosen to treat equation 1.26 for N(/i., t), as opposed to equation 
1.22 forT(/i, r), since the Q-number kernel in this case can be separated into 
a product of a function of the dummy time variable, A., times an operator which 
does not operate on the function. In this special case, the entire process of 
constructing the various local approximations can proceed leaving the second 
order differential operator untouched. The construction of the equations which 
generate the local approximations in this case is no more difficult than the 
construction presented in the previous section for the C-number kernel although 
the extraction of the solutions of these equations is considerably more difficult. 
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Since the Laplace transform of the kernel which we are considering here 
is the operator ratio of a first to a second order polynomial in p, the local 
approximation method truncates at n = 4 yielding, in that order, a partial’ differ- 
ential equation for [p. , r) which is third order in time, and which generates 
the exact solution. From equations n.37 through 11.43 we find. 


A 

11.56 

= rf I,.] (1 - V?) ^ 

“hij? 


, I 

+ 2 1^1 + (1 + 




V J 


subject to. 


0 ) 


= 0 


0 ) 


- 




0) 

'dp^ 


n.57 


On comparing, we find that we have regained equations 1.50 and 1.51 from 
paper I. which were obtained by repeated differentiation with respect to time 
of equation n.49. We have also shown in paper I. that equations 1.50 and 1.51 
are equivalent to equations 1.44 through 1.49 which were used there to obtain 
numerical evaluations of the exact quasi-linear diabatic solutions for several 
interesting initial conditions. Thus, the numerical solutions presented in paper I. 
can be viewed as the fourth local approximation to the probability distribution 
function. We further constructed the adiabatic approximation in paper I. which 
we found is governed by. 
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11.58 


dr 


= !;^l 



V 


and, in addition we presented one numerical evaluation of the solution of this 
equation for a narrow parallel beam injection. In the context of the present 
discussion, equation 11.58 and the example of a solution given in paper I. re- 
present the zero'th local approximation to the test particle problem. We have 
discussed the failure of the zero'th local approximation in paper I.. In this 
section we will present the results of a numerical evaluation of the first local 
approximation to the probability distribution function for the same initial con- 
ditions considered in paper L in order to demonstrate the improved level of 
approximation in going from n = 0 to n = 1. 

a). The Laplace modes. 

We will begin our discussion of the first local approximation applied to the 
test particle propagation problem through a discussion of the first local approxi- 
mation to the Laplace modes which we introduced in paper I.. In general we 
will show that the first local approximation to the Laplace modes and the related 
eigenvalues is a significant improvement to the zero'th local approximation for 
several of the low order Laplace modes when the Laplace variable is near zero. 

We have seen in paper I. (equation 1.64), that the Laplace modes are gen- 
erated by the eigenvalue equation, 


35 



SEPRODUCIBILITy OF'TH® 
©BiiKNiyj PAGE IS POOR 


ti ~ Cm> p) 

~ (1 - jj?') K;0, p ) + ^^(p) ypjjj., p) = 0 11.59 

'*v 

We showed generally that K{(j.^p) is a positive definite function of ^4 when p is 
real and positive (this property is explicitly demonstrated in equation n.53 for 
the special case being considered here). Thus, the Sturm-Liouville method 
could be applied to generate a complete set of non-degenerate orthogonal eigen- 
functions with a discrete set of eigenvalues. Several of these eigenfunctions 
and eigenvalues were computed numerically and presented in paper I. for various 
values of p. 

The n'th local approximation to the eigenfunctions and eigenvalues generated 
by equation n.59 is obtained by substituting the n'th convergent, (m, p), in that 
equation for (m> P) ■ The zero'th convergent to K (^t, p) is 

KqQjl, p ) = n.60 

1 -h 

This expression explicitly demonstrates a general feature of the zero'th con- 
vergent to (/li , p) which was pointed out in paper I,; i.e., the zero'th convergent 
contains a zero at = 0 which leads to doubly degenerate, discontinuous eigen- 
functions in the zero'th local, or adiabatic approximation level. The failure of 
this level of approximation was shown in paper I, to be directly attributable to 
the discontinuous nature of these eigenfunctions. 
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In Appendix B of this paper we attempt a systematic expansion in powers 
of p of the eigenfunctions and'' eigenvalues using the adiabatic quantities as the 
leading approximation and also using degenerate perturbation analysis to remove 
the degeneracy in this leading approximation. There we find unbounded first 
order corrections to the adiabatic eigenvalues due to the discontinuities in the 
adiabatic eigenfunctions. Thus, the discontinuities in the adiabatic eigenfunctions 
lead to the failure of another standard perturbation method. 

The first convergent to K p) is. 


K p) = -!■ 

1 


1 

(1 + 


P 


11.61 


This convergent is positive definite in ij, for p real and positive and therefore 
does not lead to the problems contained in the zero'th local approximation which 
were discussed above. On substituting this convergent into equation 11.59 for 
K (/i, p), we find that the low order eigenfunctions and eigenvalues are very well 
approximated for small p. In figures 6 and 7 we plot the differences between 
the exact and first local and V '2 several small values of p, and in Table I 
the exact and first local eigenvalues are compared. As in paper I., these results 
were obtained through numerical integration using a Runge-Kutta integration 
routine and the "shot-gun method" on equation 11.59 containing the first convergent 

f\j 

instead of ^ (/x, p). All of these comparisons have been made using the ex- 
ponential correlation function given by equation 11.53 
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Pigtire 6, The difference between the first Laplace mode and its first local approximation is 
shown as a function of p,. Notice that the difference has been multiplied by a factor of 
100 in this figure, and that the difference generally decreases with decreasing Laplace 
variable, p. 






Figure 7. The difference between the second Laplace mode and its first local approximation 
is shown as a function of pc. Notice that the difference has been multiplied by a factor 
of 100 in this figure, and that the difference generally decreases with decreasing Laplace 
variable, p. 



Table I 


p 

^1 

X 1 (local ) 

X 2 

X 2 (local) 

0. 1 

0.566 

0.576 

2.640 

2.675 

O.Ol 

0.2750 

0, 2753 

2.360 

2.360 

0.001 

0.1705 

0.1705 

2.317 

2.3175 


b). Reinitialiaation 


We obtain the equation which generates the first local approximation to 
N(^, r) by using the kernel given by equation 11.51 in equation 11.18. We find, 

br 

II. 6 2 

= J9qO) -rf 8 (t) (^) 

b/Ci^ b/Li^ 


b^N, 

Nj -i- j9 O) 

bu2 


in which i9^ (fi) = (1 - jj?y(n,^{jj.) and, 

HljjOLi) = Ul (^) = _ 11.63 

1 + (1 + 

We have found numerical solutions to this equation by first determining the 
reinitialized conditions from the initial conditions through, 


d^N.Oii, 0^) 

N O, 0^) + ry J9 O) i = NO, 0) 

dp? 


11.64 
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and then by using this result for the initial condition on r) to integrate 

the coupled system of equations, 


3G 


11.65 


and, 

G = N -i- 7 j2 J9 (fj.) n.66 

which are equivalent, with some further boundary-initial conditions discussed 
below, to equation 11.62 for positive r. 

Some examples of reinitialization are given in figures 8, 9 and 10 for the 
exponential correlation function. These results have been obtained by numerical- 
cally integrating equation n.64 to obtain {p, 0"^) given N(u , 0), and then dif- 
ferentiating that result with respect to p in order to obtain fi(/i, 0^). In each 
case the initial condition is considerably smoothed in the process of obtaining 
the reinitialized condition. It is our impression that this reinitialization pro- 
cedure accomplishes at least some of the smoothing which is provided by the 
initial transients in the exact solutions' of the quasi-linear diabatic equation. 

c). The probability distribution function 

The boundary-initial conditions which are necessary to integrate equation 
11.64 through 11.66 are determined as follows: From the definition of N(m, r) 
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given in equation n.50 we see that N(-l, r) = 0, and from equations 11.64 through 
11.66 we further see that N(l, r ) = N(l, 0) because (1) = (1) = 0. Since 

equation IL62 is homogeneous in N(/i, t) we choose N(l, r) = 1 without loss of 
generality. Thus, we have established the boundary conditions on N(fi , r ). 

The boundary conditions on G(ct , r ) are given by G(±l, r ) = N(±l, r ) from 
equation IL66 because (±1) = 0, and the initial condition on G(^t, r) is given 
by G(a, 0^) = N(/z, O'*') from equations 11.64 and 11.66. 

We have rewritten equations 11,65 and n.66 as difference equations, and 
have used an Implicit- explicit integration routine to determine N(^, r). Given 
N(^ , r), we compute G(^t, r + Ar) from equation H.ll, and then through an 
implicit matrix inversion, we determine N(a^, r + At) from equation n.l4. 

This procedure is repeated in order to march forward in time. We have 
used the exponential correlation ftinction in order to compute (^) 

In figure 11, we present the first of these numerical results for a parallel 
beam injection. The initial condition in this figure is a guassian in fx, centered 
at /j, = 1; it is identical to the initial condition which was used to obtain the 
quasi-llnear diabatic solution presented in figure 3 of paper I. These two 
figures should be compared carefully in order to appreciate the accuracy of the 
local approximation in this case. To facilitate this comparison, we also present 
superpositions of both solutions at one and two Larmor periods in time in 
figures 12 and 13 respectively. In each of these drawings, the continuous line 
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Figure 11. The first local approximation to the quasi-linear gjnco-phase 
averaged probability distribution function for the evolution of a narrow 
beam injected parallel to the mean magnetic field. Compare with 
figure 3 of paper I. 
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Figure 12. A direct comparison, at one Larmor period, of the quasi-linear gyro-phase 
averaged probability distribution function with its first local approximation for the 
narrow parallel beam injection of figure 11. The squares give the first local 
approximatio n . 




r 


T 


T — r 


1 — r 


1 — I — 1 — I — i — I — I — I — ! — \ — I — I — r 



Figure 13. A direct coraparison, at two Larwior periods, of the quasi-litiear gyro-phase aver 
aged probability distribution function with its first local approximation for the narrow 
parallel beam injection of figure 11. The squares give the first local approximation. 




curve represents the quasi-linear diabatic solution, and the curve made up of 
small squares is its local approjdLmation. Figure 13 does contain both solutions 
but the overlap of the solutions is so good, that the continuous line is not visible. 
Figure 13 should also be compared with figure 6 of paper I.. This comparison 
provides an example of the visual evidence which has lead us to an examination 
of the normal mode expansions of these solutions, and then, has lead us further 
to the local approximation method. It seems clear from this comparison that 
the long time structure of the probability distribution function is dominated by 
(m> P) small values of p. In this case, the effect of reinitialization is 
slight, but nevertheless, can be seen in figure 11. 

Figure 14 contains the local approximation to the probability distribution 
function for the cross field beam injection presented in figure 4 of paper I. 

In this case, the effect of reinitialization is clearly visible. Notice that, al- 
though the local approximation does not contain the dramatic initial transients 
in the quasi-linear diabatic solution, as soon as these transients decay, the two 
solutions become almost identical. A direct comparison of these two solutions 
at one, two and three Larmor periods is presented in figures 15, 16 and 17 
respectively with the same annotation as that used in the previous example. 

Once again, the accuracy of the local approximation is remarkable, and also, 
through a comparison of these figures with figure 7 of paper I. further motivation 
for the normal mode expansion can be obtained. Notice that by three Larmor 
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Figure 16, A direct comparison, at two Larmor periods, of the quasi-linear gyro-phase aver- 
aged probability distribution function with its first local approximation for the narrow 
cross field beam injection of figure 14. The squares give the first local approximation. 






Figure 17. A direct comparison, at three Larmor periods, of the quasi-linear gyro-phase aver- 
aged probability distribution function with its first local approximation for the narrow 
cross field beam injection of figure 14. The squares give the first local approximation. 




periods the probability distribution function has almost reached isotropy 
-independence) . 

In figure 18 we present the local approximation to the probability distri- 
bution function for the broad parallel beam injection presented in figure 5 of 
paper I. Once again the effect of reinitialization and the lack of initial transients 
are readily apparent. Because of the development of a numerical instability 
in the quasi-linear dlabatic solution beyond r = 1.0 Larmor period (see paper I.), 
we present only one direct comparison of these solutions at one Larmor period 
in figure 19. We have also carried the local approximation to this solution out 
to large times. In figure 20 we have plotted the local approximation every two 
gyro-periods from four to sixteen Larmor periods. As time increases, this 
solution tends toward the isotropic distribution in which the levels of the solution 
on either side of the steep transition region near = 0 become equal. There is 
'some evidence from this figure that the slope of the solution in the vicinity of 
yLi = 0 IS increasing with time. Remember that in the adiabatic or diffusion 
approximation the slope of the probability distribution function quickly develops 
an infinite value at = 0. In the next paper in this series we will show that 
this slope, both m the quasi-linear diabatic solution and its local approximation, 
approaches a large but finite value at = 0 as the time goes to infinity. It is 

the development of this large slope at = 0 which will eventually lead us to 

suspect the quasi-linear approximation to the probability distribution function 
as a viable leading approximation. 
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Figure 18. The first local approxiraation to the quasi-liaear gyro-phase 
averaged probability distribution function for the evolution of a broad 
beam injected parallel to the mean magnetic field. Compare with 
figure 5 of paper I. 
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Figure 19. A direct comparison, at one Larmor period, of the quasi-linear g 3 rro-phase aver- 
aged probability distribution function with its first local approximation for the broad 
parallel beam injection of figure 18. The squares give the first local approximation. 





Figure 20. The first local approximation to the quasi-linear gyro-phase averaged probability 
distribution function for the broad parallel beam injection of figure 18. The curves give 
the probability distribution function from four to sixteen Larmor periods in intervals of 
two Larmor periods . 




In general, it seems that we can conclude that the local approximation to 
the quasi-linear diabatic probability distribution function is very accurate after 
an elapsed time of about one Larmor period during which rapid wave-like 
transients decay away in the diabatic solution. The local approximation is 
reinitialized so that after the decay of the transients the diabatic solution and 
its local approximation join together and then evolve in unison through the 
longest times . . . several Larmor periods in some cases . . . for which our 
numerical solutions have been obtained. Of course, these general conclusions 
must be qualified by the fact that we have not even attempted a thorough in- 
vestigation of the available parameter space; we have studied a single particle 
energy for which the particle Larmor radius is equal to the correlation length 
in the turbulent magnetic field, and we have considered only a single value for 
the parameter, p , which measures the relative strengths of the random and 
mean fields. Furthermore, we have not addressed the question in this paper 
of what happens to any of these solutions for times much larger than those 
available to us in our numerical results. In the next paper in this series we will 
develop strong evidence that the local approximation remains faithftil to the 
diabatic solution for all post transient times. 

We have only mentioned briefly in paper I. the many precautions which we 
have taken to ensure the accuracy of the numerical results presented in that 
paper and also in this one. We feel that the most convincing argument that we 
can make for this accuracy is based on the direct comparison figures which 
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we have presented in this section. The completely es^licit marching routine 
which we used to obtain the quasi-linear diabatic solutions in paper I. is far 
removed from the explicit-implicit routine which we used to obtain the local 
approximations presented in this paper. If either of these integration routines 
were to develop significant errors, it would be highly unlikely that the integration 
results would compare as favorably as they do for such a variety of initial 
conditions. On this basis, our confidence in the accuracy of both integration 
routines remains high. 

V. DISCUSSION 

We have developed a new approximation method which we call the local 
approximation method. We believe that this new method is applicable to a wide 
class of equations of importance in many kinetic and transport theories. In 
the local approximation method an expansion for small values of the Laplace 
variable is carried out in order to obtain the slow evolution of the mechanical 
system under consideration. This procedure is in full accord with the aims of 
standard kinetic theory expansions which rely, however, on the existence of a 
small parameter (e = ratio of two time scales) to obtain the evolution on the 
slow time scale. The standard approach yields asymptotic expansions in the 
small parameter which often become inadequate when non-uniform behavior 
occurs in the phase space associated with the mechanical system. For the test 
particle propagation problem under consideration here a non-uniformity occurs 
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in the vicinity of yit = 0 due to mixing of the two time scales. Our new technique 
successfully removes this non-uniformity. 

We have applied the local approximation method to two examples in order 
to demonstrate its use in detail. The first of these was a simple exactly soluble 
problem all of whose local approximations were available as explicit anal 5 ^c 
expressions. The second example was the quasi-linear diabatic equation for 
test particle propagation in an axi-symmetric slab model of magnetostatic 
turbulence which was derived from first principles in paper I.. For this example 
the local approximations to the probability distribution function are available 
only numerically; we have presented only the first local approximation in this 
paper. For both examples the available local approximations generally are 
remarkably successful, but, there is one exception. 

The conditions under which the local approximations either succeed or fail 
remain essentially unknown to us. In fact, we have not formulated a quantitative 
criterion for the success or failure of any of the local approximations. For the 
two examples presented in this paper we find the local approximations generally 
successful based on visual inspections of plots of the local approximations 
compared to plots of the exact solutions. These are not asymptotic approxima- 
tions which we have constructed here, and the usual tests for asymptoticity 
cannot be applied. Often we find that the difference between the exact and 
approximate solutions is large, but in these cases the approximate solutions 
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nevertheless appear to give an accurate representation of the average behavior 
of the exact solutions. Even more often we find the approximate solutions fit 
the exact solutions extremely well. Based on the results of these two examples, 
and because this method is so unique, we feel that a much deeper investigation 
of the local approximation method is justified. 

The particular construction which we have presented in this paper is 
probably part of a more general approach. We have quasl-uniformized the 
Maclaurin series expansion of the Laplace transform of the kernel through the 
use of some of the available Fade approximants. Other of the available Fade 
approximants may be useful. We imagine also that there are other uniformi- 
zation techniques which could be used, and which would result in other approxi- 
mation methods. The exact solutions of the two examples which we have pre- 
sented here are clearly dominated for large times by contributions to their 
inverse Laplace transforms from the vicinity of the Laplace space origin. 

Thus, the approximations to the Laplace transforms of the relevant kernels 
which we have used, which are most accurate in the vicinity of the origin, are 
appropriate. In other applications, various lypes of uniformization techniques 
might be applied to Taylor series expansions about other points on the complex 
Laplace plane. It would seem to us that the key to this method is to introduce 
a simplification of the exact solution under consideration which is most accurate 
in a domain in the complex Laplace plane that has been determined, through 
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some other means, to dominate the inverse Laplace transformation. The use 
of Fade approximants for this purpose is particularly attractive since the 
result of this simplification is likely to be explicitly inverse Laplace trans- 
formable, and the failure of the simplification far away from the point of 
domination is not likely to be catastrophic. 

Based on visual inspection of the appropriate curves for our exactly soluble 
example, we would say that the first local approximation remains viable for all 
positive time. From our construction of the first local approximation to the 
low order Laplace modes, we conclude that the approximate eigenfunctions 
and eigenvalues accurately represent their exact counterparts no matter how 
small the Laplace variable. These conclusions indicate that the first local 
approximation to the quasi-linear diabatic solutions presented in paper I, may 
remain accurate beyond the ends of the available numerical evaluations. In 
the next paper in this series we investigate this issue more fully and conclude 
that this indication is correct. 
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APPENDIX A 


THE RELATIONSHIP BETWEEN EQUATION H.l AND 
OTHER KINETIC THEORIES 


Consider the pair of coupled equations, 

^ + HoF=eLG A.l 

and, 

|^+HG=IF AO 

in which, 

L = AI A.3 

where A, Hg, H, and I are linear time independent operators which operate on 
ever3d;hing to their right. The generalization in which H and I are time depend- 
ent operators will he considered later in this appendix. 

Equation A.2 can be solved formally; we find, 


G(r) = U(r) G(0) i 



dX U(>.) lF(r - X) 


A.4 


in which, 


U(r) = e"»^ 


A.5 
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Upon substituting equation A.4 into equation A.l, we find, 


^T 


+ 


H,F = 


eLU(r) G(0) + e 


dk -K(k) F(r ^k) 


A.6 


in which, 

K(\) = AI I A.7 

In the following, we will show that equation A.6 is equivalent to equation 
n.l for the purpose of appl 5 dng the local approximation method. Then we will 
demonstrate how all of the fundamental approaches to statistical mechanics 
can be put into the structure of equations A.l and A.2 with, in some cases, time 
dependent operators. Those cases with time dependent operators will be con- 
sidered last. 

In situations where the operator, Hq , appears it is the Hamiltonian or 
Liouville operator which generates the effects of an external or averaged force 
field. In paper I. equations 1.8 and I.IO contain such a term, but it does not 
appear later in that paper because of the gyro-phase averaging done in section 
III. Thus, the effects of an external or averaged force field have not been in- 
cluded in equation n.l through . 

The quantity G(0) plays the role of fluctuations or correlations in the initial 
conditions; it is traditional to ignore it. In the few cases that G(0) has not been 
set initially to zero it has been found that it has negligible effects beyond the 

7 

initial transients in the initial value solutions which are available. We feel that 
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the neglect of G (0) in situations where a non-Markovian theory is necessary 
has not been justified, however, for the purpose of this study we have conformed 
to tradition. 

It is important to realize that the inclusion of the two terms discussed in 
the previous two paragraphs does not present an obstacle for the construction 
of the local approximation method given in the main text of this paper. In the 
local approximation method the effects of the integral term in equation A.6 are 
approximated. If these two additional terms were included, they would simply 
be carried along through the calculation. Thus , for the purpose of constructing 
the local approximation method we consider equation A.6 completely equivalent 
to equation n,l, and therefore, consider equations A.l and A.2 also completely 
equivalent. 

In the following we give a list of aU the fundamental approaches which are 
basic to the construction of kinetic theory, hi each sub-section we give the 
definitions and changes of notation which are necessary to put each fundamental 
approach into the form of equations A.l and A.2. These listings are brief, but, 
the references necessary for further understanding are included. 

a). The Master equation formalism’ ’ ^ 

We consider M particles in a phase space (x ^ , v^} where i equals 1 
through M. The quantify, is the Liouville probability distribution function 
for the M particles which obeys the Liouville equation. The averaging operator, 
A, is given by. 
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A = 



d^x 

1 


The Master equation^gives the evolution of, 

F = AFM 


A. 8 


A.9 


It follows from the pair of coupled equations, equations A.l and A.2, in which, 

G = (1 - A) F“ A.IO 


H = 


M 


1 = 1 


V 

1 


1 


A.ll 


and, for particles which interact through a two-body central potential, 


M 

A<j 
j = 1 

The effects of an external force field could be included in Hq. 
b) . The BBGKY Hierarchy ’ ' ^ ° 

First consider the case in which three-body interactions are neglected. 
Then the one and two-body probability distribution functions become decoupled 
from the remainder of the hierarchy of cluster distributions. Equations A.l 
and A.2 give the> evolution of the one-particle probability distribution function, 

F ^ and the two-particle probability distribution function, F after the following 
substitutions; let the averaging operator be. 


B0 


1 j 


Bx 





A. 12 
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A.13 


A = J d^X2 

and let F = F ^ and G = F^-F^F^ Then F and G obey equations A.l and A.2 
with 


3 B 

H = V, • + V • *- I 


A. 14 


and, again for a two body central potential interaction. 


I = 


30. 


12 / 3 


3xj^ W'''i 


3 


A. 15 


In this case, even in the absence of an external force field, there is an Hq 
given by 


Ho =K = - 


3x, 


A, 16 


In the absence of a truncation of the BBGKY hierarchy as in the preceding 
paragraph, we can introduce a treatment of the entire hierarchy which is 
similar to that given in the preceding section in the Master equation formalism. 
The s -particle probability distribution function obeys the hierarchy member 
equation 


3F = 
3r 


+ K F" 

‘ S 


IF" = L F^''^ 

S S 


A. 17 
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with the operators and similar to the K, I, and L operators (L = A]) 

given in the preceding paragraph, but with summations over the appropriate 
numbers of particles. With the introduction of the M’th order vector, F, given 

by. 


F = 





p3 


A. 18 


the entire hierarchy of equations can be represented as a matrix equation, 


^ + KF - IF = LF 


A.19 


where, K,I, and L are the M'th rank matrices whose elements are all zeros 
except for the following: 


K.. = K, 


I.. = I. 


L +, = L 

1 , J +1 1 


A,20 


Now, a further averaging. A’, of equation A.19 which parallels that used to 
produce the Master equation from the Liouville equation in the preceding 
section could be formulated and would lead to a pair of coupled equations, 
equations A.l and A.2, for A' F and (1 - A')F. 

c) . The Klimentovich formalism ^ 2,13,14.15 

We illustrate the Klimontovich formalism through its application to the 
Lorentz gas. In this application the operators, H and I, which appear are time 
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independent. In a more general application this formalism contains time 
dependent operators, but then, it can be treated as in the last examples in this 
appendix. 

We consider the particle micro-density, N, for M ions and one electron, 
which obeys the Liouville equation, 


3N ^ ^ _ n 


A.21 


in which, 


(D 


ni 

0(x 


i = 1 




A.22 


and<?!> represents the two body potential interaction. The averaging operator 
in this case is given by, 


A = 



r d^v d^x 

e e 


M 

1 = 1 


A.33 


in which is the coarse grained Liouville probability distribution function 
for the M ions and the one electron. Now equations A.1 and A.2 give the evolu- 
tion of F = AN and G = (1 - A)N with 


3$ ^ 5 
3x 3v 


A.24 
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and, 


H = V - A _ I A.25 

We note in passing that this formalism is equivalent to the multiple scattering 
formalism. ^ 

d). Test particle propagation in magnetostatic turbulence with no approximations. 

This series of papers is concerned with the quasi-linear approximation to 
the ensemble averaged probability distribution function for test particle propaga- 
tion in magnetostatic turbulence. We have shown that the quasi-linear diabatic 
equation which governs the probability distribution function within this frame- 
work is in the form of equation n.l. It is also true that the equations which 
govern the exact ensemble averaged probability distribution function can be put 
in the form of the coupled pair of equations given by equations A,1 and A.2. 

Thus , the local approximation method can be applied to an exact formulation of 
the test particle propagation problem as well as the quasi-linear formulation. 

Equations 13 and 14ofKlimas and Sandri^^ can be rewritten as equations 
A.l and A.2 through the following changes in notation: Let F f , and G = f'/^’. 
Then with ensemble averaging for the A operator, and, 

Hq = K + e£ A.26 

I = - £' A.27 
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and 

H = K + eC + e' (i - A) £' A.28 

the equations due to Klimas and Sandri become identical to equations A.l and 
2 

A, 2 with {e' ) in place of e . 

In the event that the H and I operators of equation A.2 are functions of time, 
the construction of a closed equation for F is still possible, and the local approxi- 
mation method can still be applied. Here we establish the groundwork for this 
situation, and then discuss several examples in which it occurs. 

With time dependent H and I, the formal solution of equation A.2 is given 

by. 

G(r) = U(T) G(0) + U(T) r d;. U“^\) F(M A.29 

>^0 


in which. 


U(T) = T 


exp - 


ds H(s) 




A.30 


where T is the time ordering or chronological operator. The closed 
equation for F is now given by, 




+ HqF = € L U(r) G(0) + 


e 


dk K(T, k) F{k) 

•^0 


A.31 
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in. which, 


K(r, = AI(r) U(r) U"V\) I(X) A.32 

Although we have lost the-convolution form for the integral term in this equation, 
the Laplace transform can still be constructed. Rather than an algebraic equa- 
tion in the Laplace variable, p, we now obtain an integral equation in this 
variable; 


pF(p) - F(0) + HqF(p) = ^ 


f •/ — 


^ . 'V 


I(p')U(p - p') G(0) 


A.33 


2 tt 


i lO 


dco K(p, o)) F(w) 


where. 


K(p, oj) = 


2'?71 


iCO I* iCO 

dP'S(p')5^ I dp-SV) (-r^'] ^ ) 

2771 \p - pj \p +CO+P' -pj 


A.34 


'V 'V'*' 

ia which ^5i(p) and cS (p) are the Laplace transforms of, 


A(t) = I(r) U(r) I(^) 


A.35 


Of course this case is considerably more complicated than that discussed above 
in which H and I are time independent. Nevertheless, we are in a position to 
apply the local approximation method at this point by producing a continued 
fraction expansion of K(p, co) about the origin in p. 
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In the following we list a few important approaches which contain time 
dependent operators, but, by the reasoning above are amenable to the local 
approximation method. 

e) . The Klimontovich formalism^ ^ , 13 , 14,15 

When the Klimontovich formalism is applied to a more general problem 
than the Lorentz gas, then the quantity, given by equation A. 22 becomes a 
function of time through the time dependent positions of all of the particles. 
Otherwise, the notation given in section c). is left essentially unchanged by this 
generalization, and thus, equations A.l and A.2 follow. In this case, however, 
both H and I are time dependent operators. 

13, 14 

f) . The passive additive theory of Obukhov 

This is a macro-turbulence theory which gives the evolution of C, the 
concentration of a passive additive. The evolution of C is given by. 


Bt 


+ V • 


— = DV^C 
Bx 


A.36 


in which D is a diffusion coefficient, and V is the centered random field dis- 
tributed over an ensemble; it is generally a function of time. This equation 
leads to the pair of coupled equations A.l and A.2 with A an ensemble averaging 
operator, F = AC, G = (1 - A)C, 

I = - d(l-A)V|] • ^ A.37 
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H = (1 - A) V ' ^ - DV2 
ox 

A.38 

and. 


= AY • - DV^A 

OX 

A.39 

The s57mbol, [ | I ] , is used to indicate that the averaging operator does not 
operate on an 5 d;hing to the right of the symbol. 

g). Evolution of a system with a random Hamiltonian 

Equations 6 and 8 of Kaufman are identical to equations A.l and A. 2 if the 
following identifications are made: Let F = P, G = Sp, and A the averaging 
operation given by Kaufman, Furthermore, let 

Ho=<L> 

A.40 

I = - SL 

A.41 

and. 


H = Ho + (1 - A) SL 

A.42 


In this approach, SL is considered a function of time. Thus, this approach 
could be applied to the problem of particle propagation in electromagnetic 
turbulence. The approach of Klimas and Sandri was taken from this one, but 
their approach is limited to magnetostatic turbulence due to their assumption 
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of time independent operators. Kaufman's work provides the generalization 
necessary to extend the test particle analysis to time dependent electromagnetic 


turbulence. The local approximation method will still play a role in that ex 
tension of the present work. 
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APPENDIX B 


FAILURE OF DEGENERATE PERTURBATION ANALYSIS 
WITH ADIABATIC EIGENBASIS 


The equation which generates the Sturm-Liouville eigenfunctions can be 
written (see equation 1.54, paper I.) 

(D 4 ^)| m, ^>= 0 A.l 


with, 


B ^ ~ (1 - i^O, p) 

OjJ. 


2 _ 

d/j. 


B.2 


In this appendix, we adopt an indexing system for the eigenfunctions, and eigen- 
values, which differs from that used in the main text. We use the index, m, to 
denote each pair of even and odd eigenfunctions, and the index, 't, to denote which 
of each pair is being considered; the even eigenfunctions will be labeled = 0, 
and the odd eigenfunctions, = 1. We attempt a perturbation expansion of the 
eigenfunctions, and eigenvalues, in powers of the Laplace variable, p, which we 
assume is small. Thus, 

= B.3 

and 

Im. 'ty = jm, 4 p|m, 4 ••• B.4 
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The zeroth order quantities in these expansions are the adiabatic approximations. 
We have found that these adiabatic eigenfunctions are doubly degenerate with a 
pair of even and odd eigenfunctions for each eigenvalue. Thus, the zeroth order 
approximation to the eigenvalues does not depend on the index, Z . We also 
expand the differential operator, D , in powers of p according to 

D = + ... B.5 

in which 

= (-1)" ^ (1 ^ B.6 

and theJfl^ are the n'th moments of K{|U., p). On equating coefficients of increas- 
ing powers of p, we find, to zero'th order in p, 

Im, = 0 B.7 

and, to first order in p, 

|m, = - (D<1> + |m, B.8 

We expand ] m, -i using the adiabatic eigenfunctions as a basis, 

|m. = |m', ^'|m. ^*9 

and substitute this expression in to equation B.8 to obtain 

(X<«) - |m'. '“V, V\ m, 4><‘> 

B.IO 

. _ (D<i> |m, 
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From this expression, using the orthogonality of the adiabatic eigenfunctions, 
we find, 


<%, kjm, k|D<‘>|m, -ey"’] B.ll 


^(0) _>,(0) 
m n 


where 





d/Lii// 


( 0 ) 
n , k 



(1 -m") \ 




B.12 





~dfx 


(1 - \ 



by integration by parts. Notice that equation B.12 is given by zero unless “P* 
In order to prevent the divergence of equation B.ll when m = n, we set, 


k|D<^>k =0 

m , 


B.13 


and then. 




B.14 


m n 


= 0 


(m = n) 


both of which provide the first order corrections to the adiabatic eigenfunctions 
and eigenvalues. However, from equation B.12 we can see that when 't = 1, 

1 |D^i) j m, is given by a divergent expression due to the discontinuous 
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behavior of the odd adiabatic eigenftmctions atyu. == 0. Thus, in this expansion 
scheme, we find that the first order corrections to the eigenfunctions and eigen- 
values are misordered due to the discontinuous nature of the odd adiabatic 
eigenfunctions which are used as a basis for the expansion. 

The conclusions made in the previous paragraph can only be made if we 
show in addition that the degenerate eigenfunctions of defimte parity in yu. are 
the appropriate linear combination within the degenerate space (p = 0) which 
evolve continuously as p grows from p = 0. This is ordinarily done by showing 
that the off diagonal terms of the secular matrix are zero, viz., 

't 1 p D I n, = 0 m ^ 't B.15 

As argued below equation B.12 these integrals do vanish identically because 
D is even in jU. From another point of view we have graphically shown 
in paper I. (figures 6 and 7) that as p - 0 the eigenfunctions represented by 
In, of definite parity are the limiting form of the exact eigenfunctions. 
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APPENDIX C 


THE n'TH LOCAL KERNEL EOE n = ZERO THROUGH THREE 


The which appear in equation H.13 can be obtained by inverse Laplace 

transforming the convergents which appear in the boxes with dark outlines in 
the operator Pade' table shown in figure 1. From the Pade' table, 

Ko(P) = Mo C.l 


Therefore, 


KoM = MoHt) 


C.2 


Similarly, 

Ki(p) = Mg - MjP C.3 

and therefore. 


Kj(t) = Mo§(t) - Mi8(r) C.4 

At this point it appears as though we are producing truncations of a one dimen- 
sional multipole expansion of the kernel with the order of the truncation increas- 
ing with the level of the local approximation. This impression would be correct 
if we were to continue down the M = 0 column on the Pade table. However, to 
obtain the second local approximation we now leave that column and proceed 
down, the "staircase” of convergents. From this point on in our procedure we 
have found no simple characterization of the resulting K^(t). 
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The s econd convergeat is given from the Fade' table by, 




K2(p) = 


- Ml 


1 +~M-%2P 


C.5 


For large p we have. 


K2(p)P K„(oo) 


C.6 


in which. 


(CO) 




C.7 


To obtain the inverse Laplace transform of K 2 (p) we use. 


_ r\^ 


(p) = [K^Cp) -K^C®)] +K2(cc) 


_ 4(MiM 2^)^ Ml 


+ (Mo - 2 MiM"^Mi) 


p + 2M2 ^Mi 


C.8 


With this form for K„(p) find, 


K 2 (t) = (Mg - 2MiM"^Mi) S(t) + 4(MiM 2^)2 Ml e ^ l 


C.9 
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Similarly, 


K,(r) = 


Mo ^2 


S(r) 


Ml M 2 M-iM 2 j S(r) 


B,10 


+ 


27 

2 


(M2M"^)^ M 2 


g J 


We have not computed any of the higher order K^(r). 
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